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We discussed recent studies, within the framework of transport theories, on heavy ion reactions between 
charge asymmetric systems, from low up to Fermi energies. We concentrated on the analysis of ternary breakup 
events of dynamical origin occurring in semi-central reactions, where the formation of excited systems in various 
conditions of shape, excitation energy and spin is observed. At beam energies around 20A MeV, we showed 
how this fragmentation mode emerges from the combined action of surface (neck) instabilities and angular 
momentum effects, leading to the observation of three aligned massive fragments in the exit channel. At Fermi 
energies, a transition towards a prompt emission of small fragments from the neck region with larger relative 
velocity with respect to projectile and target remnants is observed. We also focus on isospin sensitive observ- 
ables, aiming at extracting information on the density dependence of the isovector part of the nuclear effective 


interaction and of the nuclear symmetry energy. 
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I. INTRODUCTION 


The behavior of nuclear matter in several conditions of 
density, temperature and N/Z asymmetry is of fundamen- 
tal importance for the understanding of many phenomena in- 
volving nuclear systems and astrophysical compact objects. 
This information can be accessed by means of heavy ion col- 
lision experiments, where transient states of nuclear matter 
spanning a large variety of regimes can be created. Actu- 
ally this study allows one to learn about the corresponding 
behavior of the nuclear effective interaction, which provides 
the nuclear Equation of State (EOS) in the equilibrium limit. 
In particular, many investigations are concentrated nowadays 
on the EOS of asymmetric matter (Asy-EOS), which has, 
comparatively to the symmetric EOS [1], few experimental 
constraints. 

We stress that the latter information is essential for the un- 
derstanding of some important properties of neutron stars, 
whose crust behaves as low-density asymmetric nuclear mat- 
ter [2, 3] and whose core may touch extreme values of den- 
sity and asymmetry. Moreover, the low-density behavior of 
the symmetry energy also affects the structure of exotic nu- 
clei and the appearance of new features involving the neutron 
skin, which are currently under intense investigation [4]. 

Over the past years, several observables which are sen- 
sitive to the nuclear EOS and testable experimentally, have 
been suggested. In this article, we will present some recen- 
t transport theory results on dissipative collisions, in a range 
of beam energies from just above the Coulomb barrier up to 
the Fermi energy regime, mainly concentrating on reaction 
mechanisms occurring at semi-central impact parameters. We 
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focus on new, exotic, break-up modes, showing how the frag- 
mentation mechanism reflects the delicate balance between 
mean-field dynamics, associated shape and volume fluctu- 
ations and rotational effects. Moreover, fragment isotopic 
properties appear quite sensitive to the low-density behavior 
of the nuclear symmetry energy. 


Il. TRANSPORT THEORIES AND EFFECTIVE 
INTERACTIONS 


Nuclear reactions are modeled by solving transport equa- 
tions based on mean field theories, with short range (2p-2h) 
correlations included via hard nucleon-nucleon elastic col- 
lisions and via stochastic forces, selfconsistently evaluated 
from the mean phase-space trajectory [5]. 

In the energy range up to a few hundred A MeV, the 
appropriate tool is the so-called Boltzmann-Langevin (BL) 
equation [6] 

EOE LIEH) =tenlfl+allf, =D 
where f(r,p,t) is the one-body distribution function, the 
semi-classical analog of the Wigner transform of the one- 
body density matrix, H (r, p, t) is the mean field Hamiltonian, 
I¢o11 is the two-body collision term incorporating the Fermi 
statistics of the particles, and 6/|f] is its fluctuating part. 

Several approximate treatments of the BL equation have 
been introduced so far. For instance, in BUU-like models, 
the fluctuating term dJ|f] is neglected [7], whereas in the 
Stochastic Mean Field (SMF) model [8] fluctuations are pro- 
jected onto the coordinate space and are implemented agitat- 
ing the spatial density profile. 

Effective interactions (associated with a given EOS) can 
be considered as an input of all transport codes and from the 
comparison with experimental data one can finally get some 
hints on nuclear matter properties. 
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Fig. 1. Three effective parameterizations of the symmetry energy: 
asystiff (dotted line), asysoft (full line) and asysuperstiff (dashed 
line). 


Unless specified, we adopt a soft isoscalar EOS (compress- 
ibility K = 200 MeV) and the free nucleon-nucleon cross 
section in the collision integral. We notice that the considered 
compressibility value is favored, e.g. from flow, monopole os- 
cillation and multifragmentation studies [9, 10]. The choice 
considered corresponds to a Skyrme-like effective interaction, 
namely SKM*, for which we take the effective mass as being 


equal to the nucleon bare mass. 


The symmetry energy Fem = Fs (kin) ov Lp ap- 


pears in the energy density functional e(p, pi) = elp) + 
p=" (p:/p)? + O(pi/p)4 + .., expressed in terms of total 
(pP = Pp + Pn) and isospin (pi = Pp — Pn) densities, from 
which the mean-field potential can be consistently derived (po 
denotes the saturation density). 

For some of the reaction mechanisms analyzed,the sen- 
sitivity of the simulation results is tested against different 
choices of the density dependence of the isovector part of 
the effective interaction. For instance, the symmetry energy 
behavior associated with three different parameterizations of 
C(p) (the asysoft,the asystiff and asysuperstiff) is displayed 
in Fig. 1, detailed description is provided in [5]. Eq.(1) is 
solved numerically by adopting the test particle method [11]. 

Within our framework, the system is described in terms of 
the one-body distribution function f, but this function may 
experience a stochastic evolution in response to the action of 
the fluctuating term ôT [f]. Then the model is suitable to treat 
the occurrence of instabilities and bifurcations of trajectories 
in nuclear dynamics. 


Ill. FRAGMENTATION AT LOW ENERGY: TERNARY 
BREAKUP 


We discuss semi-peripheral heavy ion collisions in the 
beam energy range of 10-30A MeV, focusing on the possi- 


Nucl. Sci. Tech. 26, $20509 (2015) 


197 197 
Au+ Au 
0.3 T ] T T T T T T T 


l 
b=5fm, 15A MeV (0) 
© 


0.2 


Os 


Normalized Yield 
© 


0.2 


0.1 


0 0.2 0.4 0.6 0.8 1 1.2 
Mass A (%) 


Fig. 2. (Color online) Mass distribution of the fragments Al and 
A2, and of their sum, emerging from the breakup of the PLF/TLF 
fragment, for the reaction '°’Au + °" Au at 15 and 23A MeV. Al 
denotes the smallest fragment. Masses are normalized to the average 
mass of the PLF/TLF fragment. 


ble occurrence of ternary breakup processes and on the fea- 
tures of the associated reaction products. Recent experimen- 
tal investigations have been concentrated on these exotic sep- 
aration modes, in the case of Au + Au reactions at 15 and 
23A MeV [12-14]. 

Calculations are performed with the SMF approach. Let us 
consider first the reaction at 15A MeV, b = [4—6] fm. In this 
case, we observe an almost symmetric neck rupture between 
projectile-like (PLF) and target-like (TLF) fragments, leading 
to the formation of two deformed fragments in the majority of 
the events. Important quadrupole and octuple deformations 
are noticed for the primary PLF/TLF fragments, which may 
lead to subsequent ruptures (multiple breakup) and, in partic- 
ular, to ternary breakup events, caused by surface instabilities. 
An early fragment recognition method is applied to identify 
the most probable breakup configuration of the PLF (or TLF) 
fragment. 

Figure 2 (top panel) shows the mass distribution of the 
lightest, Al, and heaviest, A2, fragments for the impact pa- 
rameter b = 5 fm. Once the masses are normalized to the 
average mass of the PLF (or TLF) fragment, the results look 
in good agreement with experimental data [12]. In particular, 
the calculations are able to reproduce the distance between 
the Al and A2 centroids and the variance of the mass dis- 
tribution. Similar results have been recently reported in the 
context of the improved QMD (ImQMD) calculations of Re- 
f. [15], where a detailed comparison with the experimental 
findings of Ref. [12] is presented. Results obtained at the 
higher bombarding energy of 23A MeV are shown in the bot- 
tom panel of the figure. The behavior observed for the mass 
distribution of fragments Al and A2 is quite close to the one 
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Fig. 3. (Color online) Calculation (BLOB) of the evolution of 


the size of inhomogeneities developing in the potential landscape 
as a function of time for the reaction *°°Xe+!*4Sn at 32A MeV 
(colours) and 25A MeV (contours) for a central impact parameter. 


obtained at 15A MeV, however variances are larger, reflecting 
the more dissipative dynamics. 


It should be noticed that, in the calculations of Ref. [15], as 
well as in SMF simulations at 15.A MeV, the lightest fragment 
(A1) emerges mainly from the neck region, thus being located 
at mid-velocity. On the other hand, in the data analysis report- 
ed in Ref. [13], the fragment with the largest parallel velocity 
(F1) has the smallest mass. In SMF simulations, this feature 
is present, but at higher energy. Indeed, at 23A MeV, owing to 
increased angular momentum effects, the PLF (or TLF) frag- 
ment rotates before it reaches its maximum deformation and 
a subsequent break-up may take place. Thus this mechanism 
allows one to test the interplay between the development of 
surface instabilities and rotational effects, which is governed 
by the properties of the mean-field interaction. As a conse- 
quence, in the case of a PLF break-up for instance, the light- 
est fragment may emerge with large positive parallel velocity. 
Thus SMF results point to the occurrence of a reaction mech- 
anism, i.e. neck rupture coupled to angular momentum effect- 
s, which could explain the experimental observation. On the 
other hand, in QMD-like calculations rotational effects could 
be missing because of the too fast reaction dynamics and the 
reduced mean-field effects [16, 17]. 


In the calculations discussed above, the sensitivity to the 
symmetry energy parametrization has not been explored so 
far. It would be extremely interesting to extend this kind of 
investigations to reactions involving neutron-rich (or even ex- 
otic) nuclei [18]. Indeed the reaction dynamics could be af- 
fected by the neutron-enrichment of the neck region, related 
to neutron skin effects and/or isospin migration mechanisms 
[5, 19]. In this case, one would also expect an important sen- 
sitivity of the reaction mechanism, and of the features of the 
emitted fragments, to the isovector terms of the nuclear poten- 
tial, opening interesting perspectives towards the extraction of 
new, independent information on the density behavior of the 
nuclear symmetry energy [20]. 
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IV. SOLUTION OF THE FULL BL EQUATION FOR 
FERMIONIC SYSTEMS 


We now move to discuss recent developments concern- 
ing the solution of the BL equation (1) in full phase space. 
This is the aim of the Boltzmann-Langevin One-Body mod- 
el (BLOB) [21]. In this approach the conventional Uehling- 
Uhlenbeck average collision integral is replaced by a simi- 
lar form where binary collisions, instead of acting only on t- 
wo test particles a, b, rather involve extended phase-space ag- 
glomerates of test particles of equal isospin A = aj, d2,..., 
B = by, bg,..., to simulate collisions between nucleon wave 
packets: 


If] + 6I[f] = of P fæ W (aB}ecDp) F(AB>cD) 


dp» do 
+ ac (Iva—wl 53) F(apscp), (2) 


where 
F(anscp) = |(1-fa)(I-fe)fofo — fafa (I-fo)(1-fo)]. 
At each time interval, the full phase space is scanned 
for collisions and test-particle agglomerates are redefined 
in phase-space cells of volume hë. The above procedure 
introduces correlations which are then exploited through a 
stochastic collision procedure. As a consequence, fluctua- 
tions develop spontaneously in phase space, with the correct 
amplitude, corresponding (at equilibrium) to a variance equal 
to f(1 — f) for volumes containing Nes test particles, Miest 
being the number of test particles per nucleon [22]. Since the 
Nicest test particles involved in collisions which are not suc- 
cessful can be sorted again in new agglomerates to attempt 
new collisions, in the same interval of time, the cross sec- 
tion contained in the transition rate W in Eq.(2) corresponds 
to the nucleon-nucleon cross section divided by Nies: © = 
onn/Nest- We notice that the transition rate W(ABscD) is 
the average of the elementary transition rates W (abs+cd) over 
the ensemble X of all the couples of test particles belonging 
to the agglomerates A and B and with velocities vg and vp. 
The BLOB model applies a precise shape-modulation tech- 
nique [23] which ensures that the occupancy distribution does 
not exceed unity in any phase-space point in the final s- 
tates; this leads to a correct Fermi statistics for the distribu- 
tion function f, in terms of mean value and variance. The 
main constraint of the above procedure is to impose a phase- 
space metric characterised by phase-space cells of volume 
h?. However, the metric of each (momentum and coordinate) 
sstochasticpace is unconstrained in genaral. In the present 
calculations we impose the maximum compactness for the 
agglomerates in momentum space which does neither violate 
Pauli blocking nor energy conservation. However, beyond the 
present application to Fermi energies, further attention should 
be paid to the compactness of the wave packets also in coordi- 
nate space when dealing with effects like collective flow and 
stopping, which become relevant at intermediate energies. 
The BLOB approach has been tested in the schematic case 
of nuclear matter contained in a box with periodic bound- 
ary conditions. For systems initialized inside the unstable 
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Fig. 4. (Color online) Multiplicity of fragments larger than Be in the reaction ‘°°Xe+'*“Sn as a function of the incident energy for a central 
impact parameter (left panel) and as a function of the impact parameter for 25, 32 and 45.A MeV (right panel). 


(spinodal) region of the nuclear matter phase diagram it ap- 
pears that the stochastic collision integral induces thermal 
fluctuations which are amplified by the unstable mean-field, 
in astochasticgreement with the standard nuclear dispersion 
relation [24]. 


V. FUSION-MULTIFRAGMENTATION COMPETITION IN 
CENTRAL COLLISIONS 


Before turning to discuss the phenomenology of semi- 
central reactions, we test the BLOB approach for the de- 
scription of central heavy-ion collisions at Fermi energies, 
which have been widely investigated in the past years, both 
from the experimental and theoretical point of view [9, 10]. 
For convenience (because experimental data exist, measured 
by the INDRA collaboration [25]), we simulate the reaction 
136Xe+124Sn at 25 and 32A MeV: in this situation the system 
is close to the threshold between fusion and multifragmenta- 
tion (below for 25A MeV, above for 32A MeV). The BLOB 
transport model is employed with the prescription of a soft 
EoS (with K = 200 MeV), a linear Asy-EoS, an in-medium 
cross section as prescribed in ref. [26], with Mies = 40, and 
with a statistics of 600 events. The dynamics is followed for 
a time of 300 fm/c. 

Figure 3 shows the evolution of the size of inhomogeneities 
which develop in the potential landscape as a function of time 
for the two incident energies 25 and 32 A MeV. These in- 
homogeneities correspond to blobs of matter which, if the 
dynamics is sufficiently explosive, may eventually leave the 
system as emitted fragments. At short times, spinodal in- 
stabilities tend to split the system in several fragments of 
comparable size, corresponding approximatly to the region 
of Neon. Later on, this mechanism enters in competition with 
a prostochasticcess of partial coalescence determined by the 
mean-field resilience, which tends to revert the system to a 
compact shape, and which is more effective for the lowest 
incident energy. 

This picture results in a large variety of exit channels 
for one given macroscopic initial condition (entrance chan- 


nel), ranging from fusion to multifragmentation and passing 
through very asymmetric configurations with a small frag- 
ment multiplicity, eventually binary, and a large size asym- 
metry. Such intermediate configurations may recall asym- 
metric fission with the difference that the distribution of mo- 
mentum transfer could reach large values, more compati- 
ble with multifragmentation. These results indicate that the 
BLOB dynamics is able to describe the development of a 
large variety of configurations in fragmentation events of 
dissipative collisions, even at the lowest energy considered 
(25A MeV), in agreement with the experimental observation- 
s. We notice that the latter result was not reproduced by SMF 
simulationstochastics, which, owing to the approximate treat- 
ment of fluctuations, underestimate fragment production in 
low-energy reactions [27]. 


VI. OVERVIEW ON A LANDSCAPE OF EXIT CHANNELS 


To present a more general survey, here we extend the set 
of simulations to a full range of impact parameters, and com- 
plete the dynamical calculation with a further decay process 
of the hot fragments through a statistical evaporation mod- 
el [28]. We still keep the same reaction, 8°Xe+1?4Sn, but 
we additionally vary the impact parameter b. 

The left panel of Fig. 4 is still restricted to central colli- 
sions and shows the evolution of the multiplicity of (prima- 
ry) fragments larger than Be as a function of the incident 
energy (coloured distribution with the average marked by a 
solid line). As expected, fragmentation is favoured as the 
beam energy increases within the energy window considered 
here. However, it should be noticed that the distribution of the 
(final) fragments, obtained after the evaporation process has 
been considered, reaches a maximum above about 35A MeV 
and a fall above 50A MeV, due to the more and more reduced 
size of the fragments surviving the evaporation stage. 

The right panel of Fig. 4 shows the dependence of the aver- 
age primary fragment multiplicity as a function of the impact 
parameter. We observe that, at large incident energies, the 
maximum corresponds to central collisions, because the mul- 
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Fig. 5. (Color online) Asymmetry of IMF (circles) and PLF-TLF 
(squares), as a function of the system initial asymmetry, for two Asy- 
EOS choices: asystiff (full lines) and asysoft (dashed lines). 
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Fig. 6. (Color online) Average N/Z of the neck fragments obtained 
in the reactions !!?:174Sn + °%:54Ni at 35A MeV, b = 6-8 fm. Two 
parametrizations of the symmetry energy have been employed: asys- 
tiff (thick black lines) and asysoft (thin red lines). Full (open) circles 
correspond to the neutron-rich(poor) system. 


tifragmentation channel dominates the reaction cross section. 
On the other hand, at smaller incident energies, the maximum 
moves to more peripheral collisions and it is reduced to events 
with three or four fragments: this indicates that around 30.A 
MeV and for semicentral impact parameters (6 fm) the sys- 
tem oscillates between binary mechanisms and ternary splits, 
where intermediate mass fragments (IMF) are formed in the 
neck region. The latter mechanism will be discussed in the 
next section. 
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VII. ISOSPIN DYNAMICS IN NECK FRAGMENTATION 


It is now quite well established that the largest part of the 
reaction cross section for dissipative collisions at Fermi ener- 
gies goes through the Neck Fragmentation channel, with IMF 
directly produced in the interacting zone in semiperipheral 
collisions, on very short time scales [29]. It is possible to 
predict interesting isospin transport effects for this fragmen- 
tation mechanism since clusters are formed still in a dilute 
asymmetric matter but always in contact with the regions of 
the projectile-like and target-like remnants almost at normal 
densities. In presence of density gradients we expect a neu- 
tron flow to the neck clusters and the effect should be larger 
for a stiffer symmetry energy around saturation [5]. This is 
shown in Fig. 5, where the asymmetry of the neck region is 
plotted for two Asy-EOS choices (soft and stiff), compared 
to the PLF-TLF asymmetry, for SMF simulations of Sn + Sn 
reactions at 50A MeV, b = 6 fm. 

This isospin migration effect can be estimated on the ba- 
sis of simple energy balance considerations. Starting from a 
residue of mass Ares and a neck of mass Ayp we assume 
that the mass A participating in the isospin exchange is ap- 
proximately equal to the mass of the neck, while it is small 
relative to the mass of the residue. This will lead to the 
asymmetry (3 + A£) of the neck, and to a total asymmetry 
bres = [B (Ares fo A) T (8 = AB)A]/Ares = B- ABA/Ares 
of the residue, with AS to be determined by minimization 
of the symmetry energy. The corresponding variation of the 
symmetry energy is equal (apart from a constant) to 


AE sym = Ares Esyn( OR) Cres F AEsyn(pr)(8 =F A8}, (3) 


where pr and pz are the densities of the residue and neck 
regions, respectively. The minimum of the variation of A Esym 
yields 


Bur z Esyn(pR) 
Bres Esyn(P1) 


(4) 


From this simple argument the ratio between the IMF and 
residue asymmetries should depend only on symmetry energy 
properties and, in particular, on the difference of the symme- 
try energy between the residue and the neck regions, as ap- 
propriate for isospin migration. It should also be larger than 
one, even more so for the asystiff than for the asysoft EOS, as 
indicated by Fig. 5. 

The isospin migration mechanism is responsible for the 
emission of rather “exotic” (neutron-rich) fragments from the 
neck region. Moreover, as discussed above, the fragment 
N/Z is nicely dependent on the parametrization employed 
for the symmetry potential. The latter effect has been recent- 
ly investigated also in the context of the BLOB model, for 
112,124Sn + 58,64Ni reactions at 35A MeV, b = 6-8 fm, for 
which recent experimental data exist [19]. Fig. 6 shows the 
average N/Z of the fragments emitted from the neck region, 
as a function of their charge, for the two reactions considered 
and two parametrizations of the symmetry energy (stiff and 
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soft). Quite neutron-rich fragments are produced, especial- 
ly in the stiff case, in good agreement with the experimental 
data [19]. 


VIII. CONCLUSION 


We have undertaken an analysis of new fragmentation 
mode, on which recent experimental investigations have been 
concentrated, which may develop in low and Fermi ener- 
gy heavy ion collisions. At low energies, the possibility of 
observing ternary breakup processes of dynamical origin is 
explored within the SMF model, looking at the concurrent 


[1] Danielewicz P, Lacey R and Lynch W G. Determination of the 
equation of state of dense matter. Science, 2002, 298: 1592— 
1596. DOI: 10.1126/science.1078070 

[2] Lattimer J M and Prakash M. Neutron star observations: Prog- 
nosis for equation of state constraints. Phys Rep, 2007, 442: 
109-165. DOI: 10.1016/j.physrep.2007.02.003 

[3] Burrello S, Colonna M and Matera F. Pairing effects on spin- 
odal decomposition of asymmetric nuclear matter. Phys Rev C, 
2014, 89: 057604. DOI: 10.1103/PhysRevC.89.057604 

[4] Piekarewicz J, Agrawal B K, Colo G, et al. Electric dipole 
polarizability and the neutron skin. Phys Rev C, 2012, 85: 
041302. DOI: 10.1103/PhysRevC.85.041302 

[5] Baran V, Colonna M, Greco V, et al. Reaction dynamics 

with exotic nuclei. Phys Rep, 2005, 410: 335-476. DOI: 

10.1016/j.physrep.2004. 12.004 

Ayik S, Gregoire C. Fluctuations of single-particle density in 

nuclear collisions. Phys Lett B, 1998, 212: 269-272. DOI: 

10.1016/0370-2693(88)9 1315-9 

[7] Li B A, Das CB , Das Gupta S, et al. Momentum dependence 

of the symmetry potential and nuclear reactions induced by 

neutron-rich nuclei at RIA. Phys Rev C, 2004, 69: 011603(R). 

DOI: 10.1103/PhysRevC.69.011603 

Colonna M, Di Toro M , Guarnera A, et al. Fluctuations and 

dynamical instabilities in heavy-ion reactions. Nucl Phys A, 

1998, 642: 449-460. DOI: 10.1016/S0375-9474(98)00542-9 

[9] Chomaz P, Colonna M and Randrup J. Nuclear spin- 
odal fragmentation. Phys Rep, 2004, 389: 263-440. DOI: 
10.1016/j.physrep.2003.09.006 

[10] Borderie B and Rivet M F. Nuclear multifragmentation 

and phase transition for hot nuclei. Prog Part Nucl Phys, 
2008, 61 (Book Series): 551-601, and refs. therein. DOI: 
10.1016/j.ppnp.2008.01.003 

[11] Guarnera A, Colonna M and Chomaz Ph. 3D stochastic mean- 

field simulations of the spinodal fragmentation of dilute nu- 
clei. Phys Lett B, 1996, 373: 267-274. DOI: 10.1016/0370- 
2693(96)00152-9 

[12] Skwira-Chalot I, Siwek-Wilczyfiska K, Wilczynski J, et al. 

Dynamical evolution of the '*’Au + '®’Au system at 15 
MeV/nucleon. Int J Mod Phys E, 2006, 15: 495-499. DOI: 
10.1142/S0218301306004429 

[13] Wilczynski J, Skwira-Chalot I, Siwek-Wilczynska K, et al. Ob- 

servation of fast collinear partitioning of the 1°” Au + 197 Au 
system into three and four fragments of comparable size. Phys 
Rev C, 2010, 81: 024605. DOI: 10.1103/PhysRevC.81.024605 


[6 


=e 


[8 


fate 


Nucl. Sci. Tech. 26, S20509 (2015) 


role of surface mean-field instabilities, dissipative and angu- 
lar momentum effects. Indeed large quadrupole and octupole 
deformation effects are developing in binary exit channels 
of semi-central reactions, which may lead to a subsequen- 
t breakup of the PLF/TLF fragments. At Fermi energies IMF 
are directly emitted from the low-density neck region. In re- 
actions between neutron-rich systems, we observe a neutron 
flow towards the neck, yielding quite exotic primary frag- 
ments. This effect, also confirmed by calculations with the 
new BLOB model, nicely depends on the stiffness of the sym- 
metry potential, allowing one to get new information on the 
low-density behavior of the nuclear symmetry energy. 


[14] Cap T, Siwek-Wilczyfiska K, Skwira-Chalot I, et al. Aligned 
ternary partitioning of the '°’Au + !°7Au system at 23A 
MeV beam energy. Phys Scripta, 2014, 89: 054005. DOI: 
10.1088/003 1-8949/89/5/054005, 

[15] Li Y, Yan S, Jiang X, et al. Dynamical effects in heavy-ion 

induced ternary breakup. Nucl Phys A, 2013, 902: 1-20. DOI: 

10.1016/j.nuclphysa.2013.02.015 

[16] Colonna M, Ono A and Rizzo J. Fragmentation paths in 

dynamical models. Phys Rev C, 2010, 82: 054613. DOI: 

10.1103/PhysRevC.82.054613 

[17] Rizzo J, Colonna M and Ono A. Comparison of multifragmen- 

tation dynamical models. Phys Rev C, 2007, 76: 024611. DOI: 

10.1103/PhysRevC.76.02461 1 

[18] Cammarata P, Colonna M, Bonasera A, et al. Sifting through 

the remnants of heavy-ion collisions for observables sensitive 

to the nuclear equation of state. Nucl Instrum Meth A, 2014, 

761: 1-6. DOI: 10.1016/j.nima.2014.05.108 

De Filippo E, Pagano A, Russotto P, et al. Correlations between 

emission timescale of fragments and isospin dynamics in '*4Sn 

+ ©4Ni and !!?Sn + °8Ni reactions at 35A MeV. Phys Rev C, 

2012, 86: 014610. DOI: 10.1103/PhysRevC.86.014610 

[20] Li B A, Ramos A, Verde G, et al. Topical issue on nucle- 
ar symmetry energy. Eur Phys J A, 2014, 50: 9-11. DOI: 
10.1140/epja/i2014-14009-x 

[21] Napolitani P and Colonna M. Bifurcations in Boltz- 
mann-Langevin one body dynamics for fermionic 
systems. Phys Lett B, 2013, 726: 382-386. DOI: 
10.1016/j.physletb.2013.08.005 

[22] Rizzo J, Chomaz Ph and Colonna M. A new approach 
to solve the Boltzmann—Langevin equation for fermion- 
ic systems. Nucl Phys A, 2008, 806: 40-64. DOI: 
10.1016/j.nuclphysa.2008.02.304 

[23] Napolitani P and Colonna M. Boltzmann-Langevin one-body 

dynamics for fermionic systems. Eur Phys J Web of Conf, 

2012, 31: 00027-00033. DOI: 10.105 1/epjconf/20123 100027 

[24] Colonna M and Chomaz Ph. Unstable infinite nuclear matter in 

stochastic mean field approach. Phys Rev C, 1994, 49: 1908— 

1916. DOI: 10.1103/PhysRevC.49.1908 

[25] Gagnon-Moisan F, Galichet E, Rivet M F, et al. New isospin 

effects in central heavy-ion collisions at Fermi energies. Phys 

Rev C, 2012, 86: 044617. DOI: 10.1103/PhysRevC.86.044617 

[26] Danielewicz P. Hadronic transport models. Acta Phys Pol B, 
2002, 33: 45-64; Coupland D D S, Lynch W G, Tsang M 
B, et al. Influence of transport variables on isospin transport 
ratios Phys Rev C, 2011, 84: 054603. DOI: 10.1103/Phys- 


[19 


paar 


$20509-6 


202306.00292v1 


chinaXiv 


EXOTIC BREAK-UP MODES IN HEAVY ION... 


RevC.84.054603 


ChinaX iv ERAT! 


Nucl. Sci. Tech. 26, $20509 (2015) 


[28] Durand D. An event generator for the study of nuclear 


[27] Bonnet E, Colonna M, Chbihi A, et al. Investigation of col- collisions in the Fermi energy domain (I). Formalism and 
lective radial expansion and stopping in heavy ion collision- first applications. Nucl Phys A, 1992, 541: 266-298. DOI: 
s at Fermi energies. Phys Rev C, 2014, 89: 034608. DOI: 10.1016/0375-9474(92)90097-4 


10.1103/PhysRevC.89.034608 


[29] De Filippo E, Pagano A, Wilczynski J, et al. Time sequence 
and time scale of intermediate mass fragment emission. Phys 
Rev C, 2005, 71: 044602. DOI: 10.1103/PhysRevC.7 1.044602 


$20509-7 


